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Abstract 

It is predicted that the fully frustrated Josephson-j unction array on 
a square lattice with alternating columns of junctions having a reduced 
critical current I^eak = /xJo(with ^ < 1 and Jo the critical current of 
all other junctions) in certain range of will under increasing of the 
external current first change from the state with zero chirality to the 
state with nonzero one at the lower critical current Id before entering 
nonsuperconducting state at the upper critical current Ic2- 
PACS numbers: 74.50.+r, 74.60. Ge, 74.60. Jg 

A square two-dimensional periodic Josephson-Junction array (JJA) in a uni- 
form transverse magnetic field, with half a flux quantum per plaquette, is a 
realization of a fully frustrated (FF) XY model Q. Its ground state is a checker- 
board pattern of plaquettes with currents flowing clockwise and anti-clockwise 
. The chirality, defined as the sum of currents around plaquette in a clockwise 
direction, has antiferromagnetic order. An interesting generalization of the FF 
XY model was proposed by Berge et al^. For the case of JJA it corresponds 
to an array in which the critical currents of alternate columns are reduced 
Iweak = /^/oiwhere Iq is the critical current of all other junctions, and fi < 1. 
See Fig.l. For /i < 1/3 the ground state has zero chirality, while for /i > 1/3 
it has a nonzero chirality. In this article I investigate the transformation of the 
ground state under an external current applied in the direction transverse to 
the weak bonds. 

When the zero current ground state for /i < 1/3 has zero chirality one 
expects that when we apply the external current the chirality will remain zero. 
The ground state with current will look like Fig.l with all 9 — 0, except there 
will be a phase difference (j> at each horizontal junction to accommodate current 
/ = /osin<^. However, for big enough current this state may become unstable 
against fluctuations leading to nonzero chirality. To check this we write down 
the energy per plaquette as a function of deviations 9i from the ground state 
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values of phases of the superconducting nodes Q : 

E = - cos(0 + 6*1 - 6*2) - cos(0 + 6*3 - 6*4) - cos(6'i - 6I3) + ^ cos(6'2 - 6*4). (1) 

After straightforward but tedious algebra we find that the nonchiral state is 
stable for currents I < Id, given by 

Id = -^Vl- 2^-3^2. (2) 
i ^ 

As we expect, Id{^- = 0) = 1, as we than effectively have an unfrustrated array, 
and Iciip = 1/3) = 0, as for /i > 1/3 even the zero current ground state is 
chiral. 

We now determine the upper critical current, i.e. the maximum current the 
array can sustain with zero voltage, but allowing for nonzero chirality, following 
the method of Bcnz et a/Q for a uniform {fi — 1) array. 

A chiral current-carrying state is shown in Fig. 2. The phase differences 
across the junctions a, /?, 7, S obey the following conditions. The sum of currents 
flowing into a node is zero: 

sin a — sin 7 — 2 sin f3 — 0. (3) 
Net current in vertical direction is zero: 

sin/3 — — fi sin S. (4) 
The sum of phase differences around plaquette is zero: 

a - 7 + /3 + (5 = 0. (5) 
The net current per junction in the horizontal direction is 



From (3)- (6) we obtain: 



/ = -(sina + sin7). (6) 

/2 ^ ^'-^"^ (7) 
— 2ax + X ' 

where a = + \), x — cos^(^^). We now maximize /, given by (7), with 
the restriction a; < 1. Straightforward, though tedious, analysis shows that for 
/X > /ic = 0.2137... the maximum current Ic2 is achieved for some value of 
j; < 1 and is given by 

1 - + 1 

h. = \^^ (8) 

For fJ. = I Eq. (8) gives the result of Benz et al^ Ic2 ~ — 1. For n < 
/Zcjthe maximum / is achieved for a; = 1, or a = 7 which is a nonchiral state. 
Correspondingly for ^ < fic, Ic2 — Id- 
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In Fig. 3 the phase diagram of the nonuniform FF JJA in the /) plane 
is shown. We see that for 0.2137... = < M < 1/3 the array as the current 
is increased wih first change from a nonchiral to a chiral state at / = Id , and 
then to a non-superconducting state at / = Ic2 > Ici ■ 

From Eqs. (2) and (7) we obtain for current / just above Id the chirahty 

X = sin a — sin 7 + sin (3 — n sin 5, (9) 

behaves as 

x~(/-/ci)i (10) 

I observed the above transition in computer simulations (see Fig. 4). For 
the simulations I used, following Falo et al |^,a resistively shunted Josephson- 
junction array, with each superconducting node having a capacitance to the 
ground. The current was injected into 16x16 array from superconducting bars 
at the edges. Free boundary conditions where used in the transverse direction. 
Very small thermal fluctuations (10~^ in units of Josephson energy h,Io/2e) are 
introduced to drive the system away from the unstable nonchiral state and to 
check the chiral state for stability. In Fig. 4(a), (b) the phases of the supercon- 
ducting nodes for the different applied currents are represented by arrows. In 
Fig. 4(a) ^ — 0.32 < 1/3, / = 0.3 < Id, and the phases in one column are iden- 
tical which is expected for zero chirality. In Fig. 4(b) the current / — 0.5 > /d, 
the phases in a column are varying which indicates a finite chirality, which is 
shown in Fig. 4(c). 

In conclusion I have obtained the phase diagram of nonuniform fully- 
frustrated Josephson array in the applied current(/)-weak bond strength(/x) 
plane. For some range of ^ it implies a transition from nonchiral to chiral 
state under increasing of the applied current. I have observed this transition in 
numerical simulations. 

I am grateful to N.Akino, J.M. Kosterlitz, and J.B. Marston for useful con- 
versations. This work was supported by National Science Foundation Grant No. 
DMR-9222812. 
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Figure 1: Nonuniform fully- frustrated Josephson-j unction array in a nonchiral 
current-carrying state. Solid lines are junctions with strong critical current, 
Iq. Dashed lines - those with weak critical current hIq- Arrows - directions of 
currents. The phases of the superconducting nodes arc indicated. In a nonchiral 
state all 9 = 0, though this state may be unstable against fluctuations of 9. 



Figure 2: Nonuniform fully- frustrated Joscphson-j unction array in a chiral state. 
Thick lines with arrows - directions of currents, a, (3, 7, and 6 - phase differences 
across junctions. 



Figure 3: Phase diagram of the nonuniform fully-frustrated Josephson-j unction 
array in the (/x, /) plane, computed using Eqs.(2) and (8). 



Figure 4: Results of numerical simulations for a 16 x 16 array with the weak 
bonds strength fi = 0.32. In (a) and (b) phases of superconducting nodes are 
represented by arrows, (a) corresponds to current / = 0.3 and is in a nonchiral 
part of the diagram (see Fig. 3). (b) / = 0.5, which is bigger than the critical 
current for the chiral transition, (c) is (b) in chirality (sum of currents around 
the plaquette) representation. 
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